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Abstract 



We investigate the evolution of parton densities at small values of the momen- 
tum fraction, x, by including resummed anomalous dimensions in the renor- 
malization group equations. The resummation takes into account the leading- 
logarithmic contributions (a^lnx)'^ given by the BFKL equation and the next- 
to-leading-logarithmic corrections from quark evolution. We present numerical 
results for the parton densities and the deep inelastic structure function F2. 
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1. Small- a? processes in QCD perturbation theory 

Perturbative QCD provides a good description of hadronic collisions at large values of 
the centre-of-mass energy and the transferred momentum pt, in the region Aqqj^ -C 
~ 5* . This description is based on the renormalization group analysis of physical cross 
sections in terms of universal parton densities and process-dependent coefficient functions. 
The parton densities fa{x, iJ?) (a = g^, qi,g, i = 1, . . . , Nf, Nf being the number of active 
flavours) are functions of the parton momentum fraction x and the hard scale /i^ ~ p?- In 
leading-twist order, they satisfy the evolution equations 

^^^^^^ = I]7a6(w,a5(/^^))/b(^,/^^) , (1) 

where uj is the moment variable conjugate to x, 

fa{uJ,fi')= f'dxX^Ux,fi^), (2) 

Jo 

the evolution kernels jabi^^, ds) are the appropriate anomalous dimensions, and the running 
coupling 05 is determined by 

Scaling violation is systematically taken into account to leading, next-to- leading, etc., ac- 
curacy in In/i^ by expanding the anomalous dimensions to fixed order in as, 

^abiuJ, as) = 7ll^ (c^) «5 + 7l? (c^) «| + 0{al) , (4) 
and using the corresponding expansions of the coefficient functions. 

At present and future high-energy colliders, a new perturbative regime opens up, char- 
acterized by small values of x (when Ag^^ -C <^ S'). In particular, plentiful data on deep 
inelastic lepton scattering at low x are coming from the HERA ep collider and fixed-target 
experiments 0]. In the small-x regime, novel dynamical effects enter as a consequence 
of the presence of the two very different large scales pf and S. The new effects are related 
to the emission of many gluons, widely separated in rapidity and disordered in transverse 
momentum. They lead to logarithmic enhancements (a^lnx)'^ in higher-loop contributions 
to the perturbative expansions of physical cross sections. Correspondingly, the anomalous 
dimensions contain singularities at — to all orders in perturbation theory, as follows 

labium, as) = E + 1:"^ + OiaUasM') (5) 

(analogous singularities show up in the coefficient functions as well). These contributions 
may spoil the convergence of the perturbative expansion at small x. Indeed, numerical 
studies in fixed-order perturbation theory have already found this instability in the x range 
accessible at HERA energies The singular terms should be summed to all orders in 
perturbation theory to obtain reliable predictions. 



Eq. (^ may be taken to define the logarithmic hierarchy at small x. We shall refer to 
the first tower of terms A^'^^ as the leading-logarithmic (L(x) ) series, the second tower B^'^^ 
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as the next-to-leading-logarithmic (NL(a;) ) series, and so on. To L(x) accuracy, small- a; 
contributions are resummed by the BFKL equation ^ . Beyond the leading logarithms, the 
study of small-x phenomena has been the object of much recent effort, including calcula- 
tional programmes to compute (9(Q;5)-corrections to the kernel of the BFKL equation @J^, 
as well as phenomenological models incorporating certain classes of sub-leading effects 
into that equation |[^-|T0[|. At this level, there is an interplay between small- a; dynamics 
and the mass singularities associated with the low transverse momentum region. High- 
energy factorization |Tl[] provides a consistent way of combining small-x resummation and 
leading-twist factorization of mass singularities. Next-to-leading results have been obtained 
in Refs. []T^,|T3[ using this method. 



The purpose of the present work is to investigate quantitatively the impact of small-x 
dynamics on the evolution of parton densities, and, in particular, on scaling violation in deep 
inelastic scattering. To do this, we use resummed evolution kernels, and make predictions 
for the parton densities by solving the resulting renormalization group equations. Note that 
higher- twist effects, such as unitarity corrections at asymptotic energies are not dealt 



with in this treatment. Also, as no attempt is made to extend the BFKL equation to sub- 
leading orders, our results are not absolute predictions, but depend upon non-perturbative 
input densities. 

The system of evolution equations (|l|) can be separated into flavour singlet and non- 
singlet components. It can be shown that the latter are not subject to small-x corrections 
of the type (|D . In what follows it will always be understood that the non-singlet equations 
are treated in the standard way using the expansion (^) up to two-loop accuracy 



The enhanced small-x contributions are all associated with the singlet sector. The singlet 
evolution equations can be written in the matrix form 

d ( fs\ ^ f Iss Isg \ I fs\ 
rfln/i2 \ fg ) V ^9S Igg I \ fg I ' 



(6) 



where the quark singlet density is defined as fs = Ynliifq, + fqi)- the next Section we 
summarize the main resummed results on the matrix kernel in Eq. (^. In Sect. 3 we match 
them with fixed-order contributions of the type (^), find the corresponding solutions of the 
evolution equations, and determine predictions for the deep inelastic structure function F2. 

In this Letter we concentrate on presenting the main results of our analysis. A more 



detailed discussion will be reported elsewhere [17 . 



2. Anomalous dimensions at small x 

At leading logarithmic (L(x)) level, QCD evolution is only fed by the gluon channel, 
and the resummation of the 0{as/uj)'' contributions to the gluon anomalous dimensions is 
accomplished by the BFKL equation . The resulting structure of the singlet anomalous 
dimension matrix is 

where the BFKL anomalous dimension jii^) is determined by the implicit equation 

^ = o:sxilL), as = CAas/TT , x(7) = 2^(1) - ^(7) - - 7) , (8) 
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ifj being the Euler ■?/'-f unction. At large values of uj/ds, the solution of Eq. (||) can be 
obtained by expanding the characteristic function xil) fo^' small values of 7. The first 
terms of the power series expansion in as for 7l(u;) are 



In the region of small x, however, smaller and smaller values oiujjas are probed, and a full 
inversion of Eq. has to be performed. This procedure leads to a many- valued function of 
the complex moment variable lo. The perturbative branch of 7l(i-o') is defined by requiring it 
to match the expansion (^) at large lo. The full solution, found numerically using Newton's 
method, is shown in Fig. 1. We see that the all-order resummation of the perturbative 
cj-poles builds up a branch-point singularity at the value loj^ = 4 In 2 0:5 ^ 2.77 0:5 in 
the complex w-plane. Correspondingly, as u decreases through ujl the BFKL anomalous 
dimension 7l((x') departs from the fixed-order behaviour: its real part rises quickly until it 
saturates at the value 7l = ^, and its imaginary part develops a discontinuity along the 
real axis. In x-space, this singularity of the anomalous dimension gives rise to a power-like 
behaviour of the gluon density at asymptotically small values of x The only other 
singularities of 7l(c<j) on the perturbative sheet are a complex conjugate pair of branch- 
points at u = (—1.41 ± 1.972)05' [0. We avoid these singularities when choosing a contour 
for the inversion of Eq. (Q) . 

At the next-to- leading logarithmic (NL(x)) level, quarks start to contribute to the 
evolution at small x on the same footing as gluons. The NL(x) contribution to the singlet 
anomalous dimension matrix can be written as 



/ Cf 

iNL = 



1nl{u^) - ^Tj 
Is 




+ 0{aUas/uj)^) (10) 



{Tf = TjiNf = ^Nf), where the corrections 7^, 7^ to the gluon anomalous dimensions are 
unknown at present, while the contributions 'jNii^) to the quark anomalous dimensions 
have recently been computed [0,^]. In the DIS factorization scheme ||18[| 'yNii^^) is given 
explicitly in terms of the BFKL anomalous dimension 7l(i-o') as[] 

.^^^(,,r.^ «5 2 + 37L-37i r^(l-7L)r^(l + 7L) p. . .... 

^^^^"'"^^ = 3-27. r(2 + 27.)r(2-27.)^^^"^ ^^^^ 

where ..(7) . | J^,^ ' ?^ ^^T! . e^P {7^(1) + T^C ^^^^^ " ^^^^ ' ^ 



r(i + 7) [-7X'(7)]j "i'"^ ^ Wo ^ x(C) j ' 

r being the Euler F-function, and x' the characteristic function in Eq. (§) and its first 
derivative, respectively. Eq. (|TT|) introduces no singularities in the tu-plane to the right of 
ujl. As in the BFKL case, at moderate values of u, > ul, an accurate representation 
of the quark anomalous dimension can be obtained by taking the first few perturbative 
contributions to Eq. ([Tl|), 

2 /^,„\3 /^,„\4~ 



^Formulae are also given in Refs.[^|3j for the MS scheme Resuhs for both schemes will be 

presented in Ref . [p7[ . 
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Note that, unlike the expansion (|]) for 7^,, all the perturbative coefficients in Eq. (0) 
are non-vanishing, suggesting an earlier departure from the lowest-order behaviour in the 
quark sector At very small values of x, the fully resummed result (|Tl]) has to be used. 

As 7l('^) approaches its saturation value at 7 = ^, the quark anomalous dimension (0) 
increases rapidly, with the behaviour 1/(1 — 27^^)^/^, leading to stronger scaling violation. 



3. Parton evolution 

We now consider the evolution as determined by the resummed kernels described in 
Sect. 2. We match the BFKL anomalous dimensions (J^) and the next-to- leading- logarithmic 
results (p!0[) with the full one- and two-loop [|l^ expressions, subtracting the resummed 



terms to avoid double counting, as follows 

7r^(-, «.) = 7i;V) as + 7l?(-) al + t( + £ i^l? . (13) 



fc=3 ^ ^ ^ k=2 ^ ^ 



We evaluate the running coupling at two- loop level from Eq. (^). Since resummed expres- 
sions for the next-to-leading gluonic contributions 75, 7,; in Eq. (p!0| ) are not yet available, 
we set them equal to their two-loop values. We call the resulting expressions the next- 
to-leading quark (NLQ(x)) approximation. For deep inelastic lepton scattering, the quark 
anomalous dimensions are more relevant, because the structure functions couple directly 
to quarks, whereas they couple to gluons via an 0(a5)-suppressed coefficient function. 

One should note that small-x resummation does not preserve the momentum sum rules 
J^alabi^^ = l,«s) = 0, because it neglects all terms other than those singular at u = 0. 
Momentum conservation can be restored to all-loop order by introducing contributions 
to the anomalous dimensions which are subdominant at — > (leaving the leading and 
next-to-leading logarithmic series unchanged) but enforce the correct behaviour at u ^ 1. 
We do this by multiplying the resummed anomalous dimensions by an overall factor of 
(1 —u). To test the sensitivity of the results to this prescription, we shall also show results 
using a 'harder' model of momentum conservation, in which we simply subtract from each 
resummed expression a constant equal to its value at u; = 1. 

We have solved the evolution equations assuming various alternative sets of parton 
distributions as input conditions at scale Qo = 4 GeV^. In this Letter we limit ourselves to 
the case of the flat MRSDO' input distributions. Flat parton distributions at low 



are suggested by the fixed-target data of the E665 Collaboration |^ , and it is of particular 
interest to see whether they could be consistent with the steep rise at small x seen at higher 
at HERA [E^ . Data at low and small x from HERA should also be available soon. 



We refer to Ref. for a more complete discussion including different inputs. 



The resulting DIS-scheme singlet densities are shown in Fig. 2 as functions of x for 
different values of the hard scale Q^. Comparison with the fixed-order prediction shows 
that the effects of the L(x) BFKL resummation set in at x-values of the order of 10~^. They 
increase as x becomes smaller, and are still moderate in the range (x^lO^^) accessible at 
present colliders. This is in qualitative agreement with earlier estimates based on models 
which seek to extend the BFKL equation to include a running coupling |]8|,[23[| . 



We see on the other hand that, in the range of x and values shown, the resummation 
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effects from quark evolution are quite sizeable, and tend to be larger than those in the gluon 
channel, in spite of their being formally sub-leading. In particular, they set in quite early in 
x, around values of x ~ 10~^. This is in agreement with expectations based on comparing 
the perturbative expansions (H) and (p^). 

It should be noted, however, that the impact of momentum-conserving corrections is 
model-dependent. The difference between the results obtained using our preferred proce- 
dure and the harder model provides an indication of the theoretical uncertainty on the 
present predictions at small x due to unknown sub-dominant contributions. 

The results described above can be used directly to obtain resummed predictions for 
the deep inelastic structure function F2. The NLQ(x) result, again for the case of the flat 
MRSDO' input distributions, is the solid curve reported in Fig. 3. The 1993 ZEUS data |^ 



are also shown for reference purposes. We see that a steep rise of F2 at small x is already 
obtained in the HERA range, starting from flat input distributions, as a consequence of 
perturbative resummation. As in Fig. 2, the next-to-leading corrections in the quark sector 
have a far more dramatic impact than the BFKL contributions. 



4. Summary 

We have performed a study of QCD dynamics at small x by incorporating resummed 
anomalous dimensions into the renormalization group evolution equations. We have used 
the resummed equations to determine predictions for the singlet parton densities and for 
the deep inelastic structure function F2. 

Our results suggest that the flxed-order evolution of parton densities which are flat at 
small X for Qq ~ 4 GeV^ is a poor approximation. Effects from leading and next-to-leading 
resummation at small x are already sizeable in the HERA region, and give rise to a steep 
increase of the structure function F2. In particular, the dominant contribution to this 
behaviour comes from the next-to-leading corrections in the quark channel. 

At the moment, the resummed predictions are still subject to large uncertainties, from 
two sources. First, the next-to-leading calculation is not complete: the gluonic contribu- 
tions beyond two-loop order are unknown. Second, subdominant contributions, such as 
those responsible for momentum conservation beyond two-loop order, are still important 
at HERA energies. 



Acknowledgments. We are most grateful to Stefano Catani for a number of valuable 
discussions on topics related to this work. We would like to thank Gavin Salam for help 
in preparing Fig. 1. RKE would like to thank the Cavendish and Rutherford- Appleton 
Laboratories for their hospitality, and the UK Particle Physics and Astronomy Research 
Council for the award of a Senior Fellowship. FH acknowledges the hospitality of the 
Theory Group at Fermilab. 



-6- 



References 

1. See, for instance, W.J. Stirling, in Proceedings of the Aachen Conference "QCD - 20 
years later", eds. P.M. Zerwas and H.A. Kastrup (World Scientific, Singapore, 1993), 
p. 387. 

2. J. Feltesse, Saclay preprint Dapnia/SPP 94-35, rapporteur talk at the 27th Interna- 
tional Conference on High Energy Physics, Glasgow, July 1994. 

3. For recent reviews see, for instance, J. Kwiecinski, Krakow Report No. 1681/PH, to 
appear in Proceedings of the Workshop "QCD94" , Montpellier, July 1994; 

S. Catani, Florence preprint DFF 207/6/94, to appear in Proceedings of "Les Ren- 
contres de Physique de la Vallee d'Aoste", La Thuile, March 1994. 

4. R.K. Ellis, Z. Kunszt and E.M. Levin, Nucl. Phys. B420 (1994) 517. 

5. L.N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338; 

E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1977) 199; 
Ya. Balitskii and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

6. V.S. Fadin and L.N. Lipatov, Nucl. Phys. B406 (1993) 259. 

7. A.R. White, Phys. Lett. 334B (1994) 87. 

8. J. Kwiecinski, Zeit. Phys. C29 (1985) 561; 

9. L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1; 
E.M. Levin and M.G. Ryskin, Phys. Rep. 189 (1990) 268. 

10. A.J. Askew, J. Kwiecinski, A.D. Martin and P.J. Sutton, Phys. Rev. D 47 (1993) 
3775, Phys. Rev. D 49 (1994) 4402. 

11. S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. 242B (1990) 97, Nucl. Phys. 
B366 (1991) 135, Phys. Lett. 307B (1993) 147. 

12. S. Catani and F. Hautmann, Phys. Lett. 315B (1993) 157. 

13. S. Catani and F. Hautmann, Nucl. Phys. B427 (1994) 475. 

14. J. Bartels, Phys. Lett. 298B (1993) 204, Zeit. Phys. C60 (1993) 471; 
A.H. Mueller, Nucl. Phys. B335 (1990) 115, preprint CU-TP-640. 

15. V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438, 675; 
G. AltareUi and G. Parisi, Nucl. Phys. B126 (1977) 298; 

Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

16. E.G. Floratos, D.A. Ross and C.T. Sachrajda, Nucl. Phys. B152 (1979) 493; 
G. Curci, W. Furmanski and R. Petronzio, Nucl. Phys. B175 (1980) 27; 

W. Furmanski and R. Petronzio, Phys. Lett. 97B (1980) 437; 

E.G. Floratos, C. Kounnas and R. Lacaze, Nucl. Phys. B192 (1981) 417. 



-7- 



17. R.K. Ellis, F. Hautmann and B.R. Webber, Cambridge preprint Cavendish-HEP- 
94/20, in preparation. 

18. G. Altarelli, R.K. Ellis and G. Martinelli, Nucl. Phys. B157 (1979) 461. 

19. W.A. Bardeen, A.J. Buras, D.W. Duke and T. Muta, Phys. Rev. D 18 (1978) 3998. 

20. A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Lett. 306B (1993) 145. 

21. A.V. Kotval, preprint Fermilab-Conf-94/345-E, to be pubhshed in Proc. Vlth Ren- 
contres de Blois, Blois, June 1994. 

22. ZEUS Collaboration, M. Derrick et al., Phys. Lett. 316B (1993) 412; 

HI Collaboration, I. Abt et al, Nucl. Phys. B407 (1993) 515, Phys. Lett. 321B (1994) 
161. 

23. E.M. Levin, G. Marchesini, M.G. Ryskin and B.R. Webber, Nucl. Phys. B357 (1991) 
167. 

24. ZEUS Collaboration, preprint DESY-94-143. 

Figure Captions 

Figure 1: The BFKL anomalous dimension 7^ as a function of the complex moment 
variable uj/as'- a) Real part, b) Imaginary part. 

Figure 2: a) Gluon and b) quark singlet density in the DIS-scheme at different values 

of (Q^ = 10, 100 GeV^) as obtained from two-loop (dashed), two-loop -|- 
L(a;) (dotdashed), two- loop -|- NLQ(a;) (solid) evolution of flat input distribu- 
tions (MRSDO') at = 4 GeV^. The dotted curves correspond to NLQ(a;) 
resummation with the harder model for momentum conservation described in 
the text. 

Figure 3: Resummed predictions for the structure function F2. The short-dashed curve 
is the one-loop fixed-order prediction. The other curves are as in Fig. 2. 
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